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Abstract. 

We present a new method of incorporating radiative transfer into Smoothed Particle Hydrodynamics 
(SPH). There have been many recent attempts at radiative transfer in SPH lITolfTTl lsl fTsll . however 
these are becoming increasingly complex, with some methods requiring the photosphere to be 
mapped (which is often of non-trivial geometric shape), and extra conditions to be applied there 
(matching atmospheres as in 101, or specifying cooling at the photosphere as in |5]). The method of 
identifying the photosphere is usually a significant addition to the total simulation runtime, and often 
requires extra free parameters, the changing of which will affect the final results. Our method is not 
affected by such concerns, as the photosphere is constructed implicitly by the algorithm without 
the need for extra free parameters. The algorithm used is a synergy of two current formalisms for 
radiative effects: a) the polytropic cooling formalism proposed by 1 12], and b) flux-limited diffusion, 
used by many authors to simulate radiation transport in the optically thick regime (e.g. |5]). We 
present several tests of this method: 

• The evolution of a 0.07 M© protoplanetary disc around a O.SM© star (0.IEIill2]) 

• The collapse of a non-rotating 1 molecular cloud (ll3l [T2ll ) 

• The thermal relaxation of temperature fluctuations in an static homogeneous sphere (iHlilil) 
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mation 
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METHODS 

The Equation of State 

SPH alone only evolves the density and internal energy of each particle: radiative 
transfer requires extra variables such as temperature and opacity. Therefore, some kind 
of prescription is required to calculate all required variables using only the density and 
internal energy. The prescription used in this work is split into two expressions: the 
equation of state which calculates temperatures and mean molecular weights, and the 
opacity law which calculates opacities. The equation of state is the same as used by 
U2il : it accounts for the vibrational and rotational energy states of hydrogen and helium, 
as well as their various dissociation and ionisation states. 



Polytropic Cooling 



This approximation uses an SPH particle's density p,-, temperature 7], and gravitational 
potential i//, to estimate a mean optical depth for the particle UM- This relies on cal- 
culating a mass averaged column density, and a mass averaged opacity. This averaging 
process allows us to account for the particle's surroundings: hot particles may be sur- 
rounded by a cooler, denser environment, increasing the effective opacity. The formalism 
guarantees optimum cooling at the photosphere - however, this formalism cannot model 
detailed energy exchange between particles. 



Flux Limited Diffusion 

This approximation handles the exchange of energy between pairs of neighbouring 
particles, using the diffusion approximation. This approximation is only valid in the 
optically thick regime, so a flux limiter is used to allow the approximation to be expanded 
to less optically thick situations. Energy flows along temperature gradients in accordance 
with the laws of thermodynamics. However, energy exchange terms are pairwise: hence 
there is no energy loss from the system. 



The Hybrid Method 

Comparing the limitations of the above two methods, it should be clear that a union of 
these two procedures should be complementary: polytropic cooling handles the impor- 
tant energy loss from the system (which flux-limited diffusion cannot), and flux-limited 
diffusion handles the detailed exchange of heat between neighbouring fluid elements 
(which polytropic cooling cannot). Indeed, polytropic cooling's inability to model the 
detailed exchange of heat between neighbouring fluid elements - and flux-limited diffu- 
sion's inability to model energy loss - allow the two methods to work together correctly, 
modelling all aspects of the system's energy budget without encroaching on each other. 



TESTS 

The Evolution of a Protoplanetary Disc 

The evolution of a O.OVA/q protoplanetary disc around a O.SMq star can be seen in 
FigureHl The density structure is consistent with the work of Mejia, Boley et al 
0]. The Toomre Q parameter is at a minimum where the cooling is at a peak, identifying 
the importance of cooling as a criterion for disc instability. 





FIGURE 1. Density evolution of the protoplanetary disc 



The Collapse of a Molecular Cloud 



The collapse of a non-rotating molecular cloud was then simulated. The spherical, 
uniform density cloud contains IMq of material (populated by 5 x 10^ SPH particles), 
and has a radius of 10^ AU (which gives a density of po = 1.41 x lO^^^gcm^^). These 
conditions were initially investigated by [3] by solving the full radiative transfer in 3D 
(with the hydrodynamics solved in ID), and were revisited by [|12|1 . As can be seen from 
Figure |2] (left panel), the initial isothermal phase (at low densities) lasts longer for the 
hybrid method, as extra cooling can occur due to the temperature gradients in the core. 



The Spiegel Test 



To study the method's time dependence, the relaxation rate of temperature fluctuations 
in a static sphere was simulated for a series of clouds of different optical depth. The lines 
in each figure show the analytic solution for the relaxation rate (0 as a function of cloud 
opacity (divided by the wavenumber, k, of the fluctuation) Ii4i|^|l2]- As can be seen, the 
hybrid method approximates the analytic solution well. 
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FIGURE 2. Left: Evolution of the central temperature as a function of central density in the collapsing 
cloud. The dotted line shows the same results for using polytropic cooling only. Right: The dispersion 
relation co for the Spiegel Test. 



CONCLUSIONS 

We have created a new radiative transfer algorithm by merging two currently existing 
methods. By doing so, the algorithm can model in detail frequency-averaged radiative 
transfer. Several tests of the algorithm against both analytic and numerical benchmarks 
have shown it to be successful in capturing the necessary physics. 
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